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AN IMPROVED METHOD FOR DESIGN OF EXPANSION-CHAMBER 
MUFFLERS WITH APPLICATION TO 
AN OPERATIONAL HELICOPTER 

By Tony L. Parrott 
Langley Research Center 

SUMMARY 

An improved method for the design of expansion- chamber mufflers is described and 
applied to the task of reducing exhaust noise generated by a helicopter. The method is an 
improvement of standard transmission-line theory in that it accounts for the effect of the 
mean exhaust-gas flow on the acoustic-transmission properties of a muffler system, 
including the termination boundary condition. The method has been computerized, and 
the computer program includes an optimization procedure that adjusts muffler component 
lengths to achieve a minimum specified desired transmission loss over a specified fre- 
quency range. A printout of the program is included together with a user-oriented 
description. 

A field test of a muffler designed with the aid of this method was conducted on a 
helicopter with a known exhaust-noise problem. When the exhaust noises of the helicopter 
with a standard exhaust system and a similar helicopter with a muffler system installed 
were compared for hover flight conditions, the muffler system was found to reduce the 
exhaust noise by approximately 11 dB(A). No significant degradation in the engine per- 
formance was observed. 

INTRODUCTION 

There is an increasing awareness of noise pollution on the part of the general public, 
especially noise caused by aircraft. In particular, users of certain types of general avia- 
tion aircraft may be compelled to restrict operations and/or modify their aircraft to com- 
ply with existing or forthcoming noise legislation that will specify upper limits on external 
noise levels (ref. 1). In addition, measurements have indicated that internal noise levels 
of 11 light, twin-engine, fixed-wing aircraft exceed the currently accepted damage-risk 
criterion for hearing loss. Analysis of this noise indicated that engine- exhaust noise was 
the primary cause of both unacceptably high cabin-noise levels and radiated far-field 
noise (ref. 2). . 



The commonly accepted solution for excessive reciprocating-engine exhaust noise 
is the installation of an appropriate reactive -muffler system. Reactive -muffler design 
efforts in the past have relied upon acoustic -transmission-line theory as a design guide. 
However, these efforts have been characterized, for the most part, by trial-and-error 
techniques. Since trial-and-error methods tend to be expensive and inflexible, there is 
a need to provide a more rational basis for muffler design. 

The aircraft-muffler designer must achieve a number of goals. The primary goal 
is to achieve a specified exhaust-noise reduction over a selected frequency range. Other 
goals include the minimization of power loss due to back pressure and the muffler volume. 
Factors of importance are the total weight, weight distribution, geometrical layout of the 
muffler and exhaust system, and the service life. Thus, to produce a viable muffler sys- 
tem requires not only a knowledge of applicable acoustic theory but also a knowledge of 
various operational constraints and their interrelationships. 

This paper describes an improved analytical design method for expansion- chamber 
mufflers and indicates how the method is applied. The method, which was developed by 
Alfredson (ref. 3), includes the effect of mean exhaust-gas flow. In the present paper the 
theory has been formulated as an extension of the standard transmission-line theory whose 
primary shortcoming has been thought to be the failure to account for mean-flow effects 
(refs. 3, 4, 5, and 6). 

The second objective of this paper is to describe a computer program that was 
originally developed by Alfredson. A discussion is also given of the application of the 
computer program to the problem of designing an expansion- chamber muffler for an air- 
craft engine. The paper concludes by presenting the results of this application to a 
helicopter. 


SYMBOLS 


Values are given both in SI Units and in U.S. Customary Units. Measurements and 
calculations were made in U.S. Customary Units. 


A k,j’ B k,j^l 

C k,j’ D k,jJ 


elements of impedance-parameter matrix 


•Dj,D 2 denominators of elements of impedance-parameter matrix 


c average sound speed in exhaust gas 

f^ ith harmonic of engine firing frequency 
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h 

I 

i 

k 

L 

M 

N 

N 1 > N 2>1 
N 3 ,N 4 J 


enthalpy 

acoustic intensity 

imaginary unit, 

wave number, 2v/X 

length of expansion-chamber component 

average Mach number 

final component of expansion -chamber muffler system 
numerators of elements of impedance-parameter matrix 


P 

Pi 

R 

S 

S i 

s 

T 

t 

V 

V i 

*1 


average total or fluctuating pressure over cross section 

fluctuating pressure at station i 

reflection factor 

cost function (eq. (29)) 

cross-sectional area at station i 

entropy 

temperature 

time 

average total velocity over duct cross section 
fluctuating velocity at station i 

transformed variable corresponding to component lengths 
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z 


specific acoustic impedance 


y ratio of specific heats 

6 irreversible pressure loss 

X wavelength 

n series multiplication 

p mass density of exhaust medium 

S series summation 

c p phase angle between incident and reflected wave 

u> circular frequency 

Subscripts: 

b branch pipe 

f presence of mean flow 

i station designation, matrix-element index 

k,j identification of impedance -parameter element 

o stagnation value 

t termination value 

Superscripts: 

(+) wave propagating in forward direction 

(-) wave propagating in negative direction 
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Notation: 


~ amplitude of fluctuating quantity 

— normalization by pc 

Abbreviations: 

sound pressure level in decibels (ref. 0.0002 dyne/cm2) 

A-weighted sound pressure level (ref. 0.0002 dyne/cm2) 
manifold pressure 
acoustic power level 
acoustic transmission loss 

MUFFLER DESIGN METHOD 

Review of Transmission-Line Theory 

An improved muffler design method has been developed which is an extension of 
transmission-line theory to include the effects of flow, termination impedance, extended 
inlet and outlet pipes, and yielding walls. The details underlying the development of this 
method are given in references 3, 4, and 5. In this section the method is described 
through the use of the concepts of transmission-line theory. Then, in the following sec- 
tion a description of the computer program and its use will be given. The reader con- 
cerned only with the application of the computer program may proceed to the section 
"Computer Program" without loss of continuity. 

The application of the electrical transmission-line analogy to acoustic filter design 
was successfully accomplished in this country by Stewart in 1923 (as discussed in ref. 7). 
The applicability of this design tool to mufflers for reciprocating engines used in aircraft 
was investigated experimentally in the early fifties by Davis, Stokes, Moore, and Stevens 
(ref. 6). Original. credit for the development of the transfer -matrix approach to four- 
terminal (or two-port) physical systems is due to Strecker and Feldtkeller (ref. 8). The' 
widespread use of the transfer-matrix approach in mechanical, fluid, and thermal sys- 
tems is largely due to the contributions of Pipes (ref. 9). In Japan this approach was 
extended to engine -exhaust muffler systems by Igarashi, Toyama, Miwa, and Arai 
(refs. 10, 11, and 12). Although transmission-line theory is well developed, it is of 


dB 

dB(A) 

MAP 

PWL 

TL 
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interest to review briefly the fundamental ideas and assumptions of this theory to place 
in better perspective the improved design method discussed in this paper. 

The key assumptions underlying transmission-line theory in its simplest applica- 
tion to muffler design are 

(1) Fluctuating pressures are small relative to the average pressure in the system. 

(2) The duct and muffler walls are rigid. 

(3) Only plane waves are propagated ^i.e., p^ = ±pcV^ ± ^. 

(4) The temperature of the medium is constant with time and spatially uniform 
throughout the muffler. 

(5) The average velocity of the medium is zero. 

(6) Viscosity effects at the duct and muffler walls are negligible. 

(7) Fluctuating quantities are decomposable into sinusoidal components of the form 

pe*" 1 . 

Operating principle . - Reactive mufflers, of which expansion chambers are a special 
case, function by causing reflected waves to be propagated back toward the acoustic 
source. In the case of expansion chambers, these reflected waves are generated by 
abrupt changes of cross-sectional area. The objective in design of expansion- chamber 
mufflers is to arrange combinations of duct area changes in the most propitious way for 
reflecting acoustic energy back into the source over the desired frequency range. The 
simplest practical expansion- chamber muffler is shown in figure 1. It consists of two 
cross-sectional-area changes separated by a length L. The amount of acoustic energy 
associated with the reflected pressure wave p^ - ) is a function of the cross-sectional 
area ratio S 2 /S 1 , whereas the length L determines the acoustic frequency at which the 
maximum reflection occurs. 

Muffler performance .- The ratio of fluctuating pressure to fluctuating volume veloc- 
ity of a forward-going plane wave averaged over a duct cross section through which the 
wave is propagating is a useful quantity for describing the noise-reduction performance 
of a reactive-muffler system and is called the analogous characteristic acoustic impedance 
of the duct, or simply the characteristic impedance when the meaning is clear from the 
context. At an isolated discontinuity in duct area, the change in characteristic impedance 
of the duct is related to the amount of acoustic energy reflected at the discontinuity . The 
total impedance change associated with a series of discontinuities may be calculated from 
an appropriate set of boundary conditions relating fluctuating quantities across the dis- 
continuities, or it may be determined from laboratory tests. When all the impedance 
changes in a muffler system are known, the input impedance at the muffler inlet may be 
calculated by starting at the tailpipe and working systematically backward to the inlet. If 
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the source impedance is known, the total radiated power of the system may be determined 
and compared with the radiated power before the muffler system is attached. This com- 
parison is accomplished by taking the ratio of the radiated power before the muffler is 
attached to that after the muffler is attached. This ratio is called the insertion loss of 
the muffler system, in analogy with the electrical filter theory. The insertion loss, as a 
function of frequency, completely characterizes a muffler system insofar as its noise - 
reducing ability is concerned. 

Source-impedance effects . - If the amplitude of the incident wave arriving at the 
muffler inlet from the source is equal to the amplitude of the incident wave when the 
muffler system is absent, then the source is said to be nonreflecting and the effective 
source impedance is equal to the duct impedance. For reciprocating-engine exhaust sys- 
tems in general, this condition seems unlikely to hold for all frequencies of interest; how- 
ever, experimental and/or analytical studies to date have not provided reliable information 
on source impedance. In its absence, the muffler-system designer is forced to rely on 
transmission loss as a measure of muffler -system performance. 

Transmission loss .- Transmission loss is the ratio of the incident acoustic power 
to the transmitted power through the muffler. It is a single number (for each frequency 
of interest) and does not permit the prediction of acoustic performance of a complete 
engine and exhaust muffler system. Furthermore, transmission loss is not directly mea- 
surable and is influenced by the termination boundary condition. Its advantages are that 
it is conceptually simple and does provide a figure of merit when comparing muffler sys- 
tems with identical terminations, hence the traditional propensity for its use. Also, to 
the extent that the source can be considered nonreflecting, the transmission loss is an 
indication of the insertion loss to be expected. 

Matrix analysis .- The algebraic complexities involved in the analysis of muffler 
performance make the use of matrix methods attractive. These methods have been taken 
over from transmission-line theory and can be easily adapted to high-speed computers. 
Also, the concept of the "impedance-parameter matrix” provides a better understanding 
of system component interaction. The impedance-parameter matrix relates the total 
fluctuating pressure and velocity at one station in a muffler system to the total fluctuating 
pressure and velocity at a different station in the system. For example, the inlet and out- 
let total fluctuating pressures and velocities for the expansion chamber of figure 1 (i.e., 
the sum of the respective forward- and backward-going wave pressures) can be related 
by a matrix equation of the form 


— J. 





Pk" 


" A k,j 

B k,j 

V 

_ v k 


_ C k,j 

D k,[ 
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The elements of the square matrix can be related to the impedances of the system at the 
stations where the fluctuating pressures and velocities are observed. Such a matrix is 
built up from similar matrices characterizing the behavior of the muffler elements. For 
instance, continuity of pressure and volume velocity for the simple area discontinuity 
shown in figure 2(a) gives 


02 


1 

0 

Pi 




S 1 


v 2 


0 

S 2 _ 

Vl 


( 2 ) 


and for a section of constant-diameter pipe of length L, the appropriate relationship is 


02 

Vo 


cos kL jpc sin kL 


— sin kL 
pc 


cos kL 


Pi 

Vi 


(3) 


For a branch component (see fig. 2(b)) that can be treated as a lumped branch impedance 
Z b , the fluctuating pressures and velocities are related by 


(4), 


03 


1 0 

Pi 


— 

_1_^2 




pc So Si 


v 3 

_ _l 


w X 

Jt s ! 

Vl 


If there are N sections in series, then the fluctuating pressure and velocity at the inlet 
section N can be related to the pressure and velocity at the outlet section 1 by an 
impedance-parameter matrix which is equal to the product of all the component matrices 
as follows: 


A N,1 

b n,i 

N-l 

= n 

Vi,i 

B i+l,i 

5n,1 

d n,i_ 

i=l 

^i+l,i 

B i+l,i 


Transmission-loss calculation .- For transmission-loss calculations it is convenient 
to express the fluctuating pressure and velocities throughout the system in terms of the 
fluctuating pressure amplitude associated with the incident wave at the tailpipe termina- 
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tion By using the boundary condition for the radiation impedance at the tailpipe 

exit 

Pj = pcZ t Vj (6) 

the fluctuating pressure and velocity at the exit may be expressed as 



The radiation impedance Z t may be calculated (ref. 13) or measured; however, assum- 
ing Zp is known, equation (7), together with the impedance -parameter matrix given by 
equation (5), implies a knowledge of the fluctuating pressure and velocity at the muffler 
system inlet from which the pressure associated with the incident wave at the muffler 
inlet can be shown to be simpily 

P N^ = |( Pn + PCVn ) ^ 

' S 

The transmission loss is given in terms of the areas and incident-wave pressures at the 
muffler system inlet and termination as 



Since P^ + ^ is proportional to P^ + ^, the transmission loss is independent of Pj^ 
and any convenient numerical value, say unity, may be assumed for p^ + ^ for computa- 
tional purposes. 

Inadequacies of transmission-line theory .- The laboratory data of reference 6 indi- 
cated good agreement between measured and calculated transmission losses for expansion 
chambers applicable to aircraft engines. These data were taken at room temperature, 
with no mean flow, and with a loudspeaker capable of producing a pure-tone sound pres- 
sure level of approximately 140 dB in the muffler inlets. Also, nonreflecting terminations 
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were used for the muffler outlets. When similar mufflers were tested on an aircraft 
engine, the exhaust-noise reduction fell short of that predicted from the transmission- 
loss calculations (ref. 6). The authors of reference 6 believed that the discrepancies 
were related to mean-flow effects and the violation of the small-amplitude assumption of 
linear acoustics. The recent work reported in references 3, 4, 5, and 14 suggested that 
for expansion- chamber mufflers, the nonlinear effect due to large acoustic pressures was 
negligible for design purposes. Further investigations reported in references 3, 4, and 5 
suggested that the effect of mean gas flow on the termination boundary condition and other 
impedance discontinuities was the most significant factor neglected in the standard 
transmission-line theory. Consequently, mean-flow effects were included in the analysis 
of expansion- chamber performance. The results of this analysis will be summarized in 
the following section. 

Effect of Mean Gas Flow on Muffler Performance 

The improved muffler design method described herein involves the effect of the 
mean flow on the impedance changes at duct-area and branch discontinuities. Also, the 
effect on the tailpipe termination impedance is included. The effect of flow on the inter- 
nal impedance changes is a straightforward extension of transmission-line theory using 
the linearized equations for energy, mass, and in some cases momentum conservation 
across area discontinuities. The single most significant effect of the flow is associated 
with the radiation characteristics of the tailpipe. 

Tailpipe radiation . - In a duct with a termination characterized by the unflanged-pipe 
reflection factor R, the relation between the incident-wave pressure amplitude p^ and 
the reflected-wave pressure amplitude p^~^ is given by 

p(-) = Re^pW (10) 

where <fi is the phase angle between the incident and reflected wave. Thus the net 
intensity in the pipe is given by 



For a mean flow with 


(12) 


where p^ is the pressure amplitude of the forward-going wave. 
Mach number M, the acoustic intensity becomes 




2pc d 1+M)2 


(1 - M) 2 R 2 (M) 
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In this equation (see ref. 3) the reflection factor becomes a function of Mach number. An 
increase in acoustic radiation occurs as a result of the mean flow. To appreciate this 
effect, consider the ratio of If to I expressed in decibels: 


increase in PWL = 10 log 5 = 10 log (1 - M) 2 R 2 (M) 

1 . i - r 2 (0) 


(13) 


where R(M) is now a reflection factor which depends upon the Mach number of the mean 
flow. The measurements of reference 12 suggest a 3- to 5 -percent increase in the 
reflection factor and no change in phase angle for Mach number and wave-number ranges 
typical of exhaust systems. Thus R(M) in the numerator of equation (13) can be 
replaced by the reflection factor for zero mean flow R(0) with a maximum error of 
approximately 2 dB. Thus, to a good approximation, the increase in the radiated acoustic 
power level relative to that for zero mean flow is given by 


(1 + M) 2 - (1 - M) 2 R 2 (0) 

Increase in PWL 85 10 log (14) 

1 - R 2 (0) 

This equation is plotted in figure 3 with the Mach number as a parameter. From this 
plot it is seen that the radiated power increases significantly relative to that for zero 
mean flow when the reflection factor approaches 1 and as the mean-flow Mach number 
increases. Typical exhaust-pipe diameters, exhaust-gas sound speeds, and exhaust- 
noise frequencies are such that the reflection factor is greater than about 0.975. For a 
typical exhaust-flow Mach number of 0.1, the plot of figure 3 implies an increase of 8 dB 
or more in the radiated acoustic power over that for zero mean flow. 

Area contraction .- References 3 and 5 indicated that a one-dimensional description 
of wave propagation on either side of a duct discontinuity would be adequate for muffler 
design provided that the mean-flow effect was included. For flow area contractions the 
isentropic energy-balance equation was linearized to obtain one equation relating fluctu- 
ating pressure and velocity on either side of the discontinuity. This resulting equation is 
given as 


p 2 + PCM 2 V 2 = Pj + pcMjV^ 


(15) 


where the subscripts 1 and 2 denote the downstream and upstream stations, respec- 
tively, as indicated in the schematic of figure 4(a). From mass conservation across the 
area discontinuity a second linearized equation is found to be 
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Sg^^2P2 PcS^Vg — PcS^V^ 


(16) 


If these equations are formulated in a matrix format and solved for the upstream fluctu- 
ating pressure and velocity, there results 


(17) 


Thus, the effect of the mean flow is to introduce a new element into the impedance- 
parameter matrix and a correction for the element which relates the velocities. Clearly, 
when the Mach number becomes zero, the matrix reduces to the form of that in equa- 
tion (2). 

Area expansion . - Experimental evidence discussed in reference 3 suggests that the 
energy relation connecting fluid state properties across a sudden flow area expansion is 
nonisentropic. However, an adiabatic condition can be assumed to hold. Hence the stag- 
nation enthalpy would be conserved across such a discontinuity; that is, 



dh„ = T ds + — + V dV = 0 
° P 

The entropy term T ds has been shown in reference 15 to be given by 


(18) 


T ds 


s d5 
p(y - 1) 


(19) 


where 6 can be regarded as an irreversible fluctuating pressure loss due to the non- 
isentropic nature of the process. Also, the isentropic relation between the fluctuating 
density and fluctuating pressure no longer holds, but instead is replaced by 
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(20) 


By using the relation (19), equation (18) can be integrated across the discontinuity from 
the downstream to the upstream side, as indicated in the schematic of figure 4(b). There 
results after linearization 


p 2 + PcM 2 V 2 = Pl + pcM 1 V 1 - — 


(21) 


From mass conservation and the use of equation (20) to eliminate the fluctuating density, 


S2M2P 2 + pcS 2 V 2 = SjMjPj + pcSjVi + S^M^6 


(22) 


Because of the introduction of the unknown quantity 6, a third equation must be derived. 
This can be obtained from momentum conservation across the discontinuity, which gives 
after linearizing and the use of equation (20) 

+ S 2 M 2 2 jp 2 + 2pcS 2 M 2 V 2 = + 2pcS 1 M 1 V 1 + S 1 M 1 2 6 (23) 

where the momentum source term due to the discontinuity has been taken to be 

p 2 (Si - S 2 ), which is based on experimental evidence indicated in reference 3. When this 

equation is solved for 6 and substituted into equations (21) and (22), the upstream fluc- 
tuating pressure and velocity can again be related to the corresponding downstream 
quantities by an impedance-parameter matrix whose elements are given by 


1 + 


— (y - 1) - 2y 


+ y 


Mi 


k ,j 


Dl 


(24a) 


B 


k,j 






Mi 


Dl 


(24b) 
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c 


k,j 


_i_ 

pc 


y !i!i 

r S 2 \ S 2 


Di 



(24c) 




D 


S 2 s 2\ s 2 


k,j 




Di 


Mi 


(24d) 


where 


Dj = 1 + 


(y - l) 

i 


\ S 2/ 



2y 


“N 

\ 

/ \2/ \ 

vjfMi 2 + 

/SfWsi \ 

\®2/ \ S 2 J 


M, 


(24e) 


Again, these equations are seen to reduce to a form equivalent to equation (2) when the 
Mach number approaches zero. 


Area contraction with a branch .- Practical expansion-chamber mufflers usually 
have extended inlet and outlet pipes, as shown in the schematic of figures 4(c) and (d). 

Such geometrical arrangements perform acoustically as a side branch with impedance 
Zjj. The side branch is indicated as region (¥) in figures 4(c) and (d). The value of Z^, 

which is specified at the junction of regions (^T), (T), and (jl), is controlled by the rigidity 
of the end wall and by the branch length. The branch impedance Z ^ is related to the 
reflection factor and phase angle 0^ at the branch entrance by the equation 



1 +R b e 
1 - R b e 


3^b 

^b 


(25) 


If the reflection factor and phase angle at the end cap are denoted by R and 0, respec- 
tively, then R^ = R for no energy dissipation between the entrance and end cap. Also, 
0^ = 0 - 2kL, where L is the branch length and the coordinate direction is that shown 
in figures 4(c) and (d). By the appropriate selection of the various branch lengths in a 
muffler system, the transmission-loss characteristic can be tailored to specifications 
within certain limits. For a branch impedance associated with a flow area contraction 
shown schematically in figure 4(c), the impedance-parameter matrix for the simple area 
contraction can be generalized to include the branch impedance. The result is 
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1 S 1 S 2 


1 _ _±__±_“ Mi - M 


So S 


P3 


3 °3 


1 - — Mi 

S3 1 



pcs 


1 - -4 



1 S-i So C 

Mi — — Mi' 
1 -rr o_ c 1 


. . (fif 

\ S 3/ 


Z. s 3 s 3 


Mi' 



1 S 2 1 



s 3 U b s 3 ; 1 




v i 


( 26 ) 


Note that if Z b approaches infinity (i.e., the end wall becomes rigid and the branch 
length becomes zero), then the impedance-parameter matrix reduces to that of equa- 
tion (17). Furthermore, if the mean-flow Mach number approaches zero, then the result 
becomes equivalent to that of equation (4). 

Area expansion with a branch . - If the muffler inlet pipe is extended, then a branch 
impedance becomes associated with an area expansion. When this impedance is included 
in the boundary condition for an area expansion, the elements for the impedance- 
parameter matrix relating upstream and downstream fluctuating pressures and velocities 
become 


m 1 S 1 S 2 
N 1 + ~ 

Z h b 3 b 3 


k,j 


Mj 2 + (y - ljMj 3 +gi M x 4 

3 

Di + D2 


(27a) 




B k,r 


Di + D2 


N 1 S 2 L __ 2\ 

N 3 — ^-1 + yM x ) 

pcZv, S 3\ ' 


C k,j = 


PcZ b °3 

Di + D2 


(27b) 


(27c) 


15 



(27d) 



N, 


D 1 + E>2 


where Nj, N 2 , N 3 , and N 4 are the respective numerators of equations (24a, b, c, 
and d), and Dj is given by equation (24e). The quantity D 2 is associated with the 
branch impedance Z b and is given by 


»2 


_l_f 2 /fl 

Z b S 3 \S 3 



(27e) 


Again, when Z b becomes large, these equations are seen to reduce to equations (24), 
as expected. It was found by the authors of references 4, 5, and 6 that yielding walls 
could significantly affect muffler performance; consequently, in references 4 and 5 the 
impedance of the end caps was altered in a trial-and- error procedure to force the mathe- 
matical models for the branched elements to conform with measurements. For a flow 
area contraction the reflection factor and phase angle of the branch impedance at the end 
cap were taken as 0.8 and 0.01 radian, respectively, whereas for flow area expansion the 
reflection factor was taken to be 1.0 and the phase angle as 0.01 radian. The lower 
reflection factor for area contraction was used because the isentropic flow condition 
across that type of discontinuity does not allow the occurrence of acoustic energy dissipa- 
tion, which in fact is known to exist but to a much lesser extent than for an expansion type 
of discontinuity. This completes the formulations of the impedance-parameter matrix 
for the various discontinuities encountered in a muffler system consisting of a series of 
cascaded expansion chambers. 

Finite-length section .- To relate the fluctuating pressure and velocity at the 
upstream side of a discontinuity to the fluctuating pressure and velocity at the down- 
stream side of the next discontinuity joined by a constant-diameter duct (see fig. 5) car- 
rying a mean flow with Mach number M 3 , equation (3) is modified as follows: 


- 


p 4 


V 4 

/ 


L_ 


cos 


kL 


\i - m 3 ‘ 


3 kL 


sin 


jpc sin 


kL 


\1 - m 3 ‘ 


cos 


kL 


\1 - M37 \1 - M 3 ‘ 


(28) 
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Clearly, for the typical mean-flow Mach numbers of 0.1 encountered in exhaust systems, 
this correction is relatively insignificant. 

Model comparisons .- Figure 6 shows a comparison of transmission loss calculated 
by the theoretical model described in this paper with the standard transmission-line 
theory discussed in reference 6. The muffler configuration chosen for illustrative pur- 
poses is shown in the sketch in the figure along with relevant geometrical data. Note that 
there is a significant reduction of transmission loss in the frequency range 200 to 350 Hz 
due to both the mean-flow effects and the yielding end cap (reflection factor 0.8 for flow 
area contraction). On the other hand there is an increase in transmission loss above 
350 Hz over that for no mean flow. 

COMPUTER PROGRAM 

A computer program called EXRSIL has been written that incorporates the effects 
of the mean flow as well as estimated reflection factors for the muffler end walls associ- 
ated with flow area contractions. This program will enable a user to design expansion- 
chamber mufflers of up to four stages. To increase the usefulness of the program as a 
design tool, an optimization subroutine has been included that will adjust all muffler com- 
ponent lengths to approach a specified minimum transmission-loss characteristic within 
imposed component-length constraints. The physical meanings of the inputs are 
described in this section and the detailed formats are described in appendix A. Appen- 
dix B contains a listing of the program. 

Program Inputs 

Input quantities for the computer program consist of the following: 

(1) Desired minimum transmission loss at three wavelengths 

(2) Acoustic sound speed and wavelength in the exhaust system 

(3) Gas-flow Mach number in exhaust tailpipe 

(4) Optimization attempts 

(5) Number of expansion chambers 

(6) Length constraints and initial lengths for all expansion- chamber components 

(7) Cross-sectional areas of expansion chamber, annular region, and tailpipe 

Desired transmission loss .- The desired minimum transmission-loss character- 
istic in decibels is deduced with the help of narrow-band analysis of far-field exhaust 
noise over the relevant operating range of engine load and rotational speed. The engine 
should preferably be equipped with an exhaust pipe of the same length as the final com- 
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plete muffler system. A representative narrow-band spectrum of exhaust-noise- 
dominated helicopter noise is shown in figure 7. Prominent engine-bank firing frequen- 
cies are labeled in the figure as f^, f£ , and so forth. The test conditions for this spec- 
trum will be discussed in the section "Application to Helicopter." The frequency range 
over which the muffler is to be effective and the amount of transmission loss desired at 
each frequency are largely a matter of engineering judgment. Since the low-frequency 
noise components are less objectionable, it is important to transform the narrow-band 
noise levels to A-weighted levels to reduce the required minimum transmission losses at 
the lower frequencies. To see how the minimum desired transmission loss is specified, 
the bar graph of figure 8 should be studied. This bar graph represents the important 
exhaust noise components of the spectrum shown in figure 7. 

In figure 8, the shaded bars represent the abstracted narrow-band levels associated 
with the engine firing frequencies shown in the measurement represented in figure 7. The 
overall level associated with the exhaust noise was found to be 96 dB. The solid bars 
represent the A-weighted levels for which the overall level was found to be 82 dB(A). The 
arrows represent the desired levels to which the lowest, intermediate, and highest fre- 
quency A-weighted components are to be reduced. These desired exhaust noise levels are 
decided upon from a determination of the levels of other dominating noise sources and, 
possibly, other considerations such as volume and weight. The computer program con- 
structs by linear interpolation the desired transmission-loss characteristic as a function 
of log 10 X, as shown by the solid line in figure 9. This interpolation is based upon three 
pairs of numbers consisting of transmission losses and wavelengths corresponding to the 
firing frequencies indicated by the arrows in figure 8. The dashed line in figure 9 is the 
computed transmission loss for an optimized muffler which will be discussed in the sec- 
tion "Application to Helicopter." 

Wavelength and Mach number .- The wavelengths at which the designer specifies a 
desired transmission loss must be calculated from the speed of sound in the exhaust gas 
estimated from the mean temperature distribution throughout the proposed muffler sys- 
tem. Gas-flow Mach number must likewise be measured or estimated for the tailpipe of 
the proposed muffler system. 

Optimization attempts . - This input specifies the number of different directions of 
the gradient vector taken in cost-function space (see section "Optimization Subroutine") 
to achieve a minimum cost function. If a zero value of the cost function is achieved, then 
the number of optimization attempts should be increased and/or adjustments to geometric 
constraints considered. 

Geometrical inputs . - The remainder of the inputs concern the geometry of the muf- 
fler system. The volume and shape constraints are specified by the number of expansion 
chambers, the muffler stage component lengths, and the cross-sectional areas of the tail- 
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pipe, chamber, and annulus. It should be noted that the cross-sectional areas cannot be 
varied within a stage or between stages. There are four component lengths per muffler 
stage as indicated in figure 10. Constraints on these lengths are specified independently 
for each stage together with an initial length for each stage component. 

Program Outputs 

Program outputs are the following: 

(1) Input data listing 

(2) Computer-generated desired minimum transmission-loss characteristic 

(3) Record of optimization attempts 

(4) Optimized component lengths 

(5) Computed transmission-loss characteristic 

Input data listing . - This section of output simply tabulates the minimum, initial, and 
maximum values of all component lengths working from the tailpipe to the inlet. Also 
given are the cross-sectional areas and tailpipe Mach number. 

Desired transmission loss .- The minimum desired transmission-loss character- 
istic is listed along with the frequency and wavelength as determined from the mean tem- 
perature of the exhaust gas. 

Optimization attempts .- The next listing shows a record of the optimization attempts 
which consists of the cost functions and the corresponding component lengths at a partic- 
ular point in the optimization process; The message NEW POINT indicates the change of 
direction of the gradient vector in the cost-function space. 

Optimized component lengths .- This listing shows the component lengths for the 
final configuration together with the final cost function. 

Computed transmission loss ..- This section of output lists the computed transmis- 
sion loss of the final configuration at the same wavelengths and frequencies as was done 
for the desired transmission loss. 

Engineering Judgments 

Reduction of low-frequency exhaust noise generally deteriorates with decreasing 
expansion-chamber length. Thus it is mandatory to establish the maximum muffler length 
available if low-frequency noise is a problem. Most practical muffler systems will incor- 
porate flow-reversing bends to achieve necessary compactness, as shown in the schematic 
of figure 10. The transmission loss of such bends is estimated to be less than 2 dB when 
transmitting into a nonreflecting termination over the frequency range of interest in 
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reactive -muffler design. Thus, the effect of smooth bends should be relatively unimpor- 
tant for design purposes. 

The muffler geometry also affects the back pressure at the engine exhaust ports. 
Since an increase in mean back pressure has been found to correlate roughly with degra- 
dation of multicylinder -engine performance, mean back pressure is conventionally used 
as an indicator of possible loss in engine performance. Mean back pressure is minimized 
in purely expansion-chamber mufflers by making bends as smooth as possible and by 
alining inlet and outlet pipes, as shown in figure 10. No attempt is made in this program 
to compute magnitude of mean back pressure. 

All component-length constraints should be mutually consistent and compatible with 
the space requirements. If some component lengths must be held constant, constraints 
must still be specified to bracket the initial length so that a significant variation cannot 
occur. Maximum lengths of internal pipes should not be such that their separation is less 
than 1 pipe diameter and preferably 1 chamber diameter. 


Optimization Subroutine 

The transmission loss of a muffler system is treated by the optimization subroutine 
as a function of 4N variables, where N is the number of muffler stages. By systemati- 
cally changing the various component lengths and comparing the resulting transmission 
losses with the desired minimum transmission losses at 151 equally spaced values of the 
logarithm of the wavelength, the subroutine attempts to minimize a cost function. This 
cost function is defined by 


S ( L 1’ L 2’- 



(29) 


where (TL^Aj) is the desired transmission loss at the wavelength A i and (TL) c (Aj) 
is the computed transmission loss at the same wavelength for some given configuration. 
If (TL)^Ai] = (TL) c (\i) for a particular wavelength, then that term is omitted. 

The cost-function independent variables 1^ must all have boxlike constraints 
imposed on them, that is, 


L min < L i < L max 


(30) 


On the other hand, the optimization procedure must operate with independent variables 
that are unconstrained. This condition is fulfilled by defining a new set of variables 
by means of the transformation 
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siAj = ill L ™” (31) 

Lmax _ L m i n 

The details of the process by which the procedure systematically alters the com- 
ponent lengths are discussed in reference 3 and will not be given here; however, the 
process is basically a modification and extension of the steepest-descent method that 
results in an improved rate of convergence. From the user's standpoint, he can regard 
the subroutine as a minimization process that operates by following a prescribed direc- 
tion vector in the cost-function space until a minimum is achieved; then the direction 
vector is changed or updated whereupon a search for a new minimum will begin. This 
event is signified in the output listing by the heading NEW POINT. The number of new 
points (or optimization attempts) to be taken is set by the user by reading in a value of 
the integer variable NO, which is usually taken to be between 5 and 10. It is the nature 
of the optimization procedure to find the minimum cost function which is nearest the ini- 
tial configuration of component lengths. Since other more significant or preferable min- 
imums may exist, several different sets of initial configurations should be tried. At the 
end of the last minimum search, a normal return of control to the main program will be 
indicated by the message SEARCH COMPLETED, NORMAL RETURN. 

This completes the description of the computer program and its use. In the 
remainder of this paper an evaluation of the results obtained upon applying the program 
to the design of a muffler for a helicopter will be given. 

APPLICATION TO HELICOPTER 

The purpose of this section is to describe the design and to present the measured 
noise -reduction performance of a muffler for a helicopter engine exhaust. The muffler 
design was accomplished with the aid of the previously described computer program. 

Apparatus and Methods 

Test helicopter .- A photograph of the test helicopter for which a muffler system 
was designed is shown in figure 11. The physical characteristics are given in table I. 
The aircraft was powered by a six-cylinder, horizontally opposed engine. The exhaust 
gases from each bank of three cylinders are routed through a manifold and then out a 
straight stack to the side and rear of the helicopter, as shown in the photograph of fig- 

N c x rot. speed 

ure 12. The engine fundamental firing frequency is defined to be — for a 

120 

four-stroke operating cycle, where N c is the number of cylinders. For the manifold 
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geometry of this particular engine, it was convenient to correlate the exhaust noise com- 
ponents with harmonics of the engine -bank fundamental firing frequency, in which case 
N c was taken as 3. 

Acoustic measurements .- Initial base-line acoustic measurements were conducted 
for the helicopter with the engine equipped with the standard exhaust manifold system 
shown by the photograph of figure 12 and hereinafter called the standard helicopter. The 
acoustic signals emitted by the helicopter were measured with commercially available, 
piezoelectric, ceramic-type microphones and were recorded by a multichannel, frequency- 
modulated tape recorder at a tape speed of 0.76 m/sec (30 in/sec) and a center frequency 
of 54 kHz. The frequency response of the complete measurement system was within 3 dB 
over the frequency range 12 Hz to 12 kHz. 

The microphones were positioned at a radius of 30.5 m (100 ft) from the aircraft 
every 18° over one quadrant. A complete 360° noise survey of the helicopter noise was 
obtained by hovering the helicopter first at 0°, then 90°, 180°, and finally at 270°. The 
measurements were conducted over flat terrain at the approach end of runway 17 at 
Langley Air Force Base. 


Muffler Design and Installation 

Design procedure . - The computerized analytical muffler design method and required 
engineering judgments have already been described elsewhere in this report. In this sec- 
tion specific design considerations for the test helicopter will be discussed. 

Weight and volume limitations indicated that a single muffler that would handle the 
exhausts from both engine banks would be preferable to a muffler for each bank. Also, 
it was anticipated that the Y-connector used for combining the two streams of engine 
exhaust gas would provide some noise reduction of the engine odd-harmonic firing fre- 
quencies and some amplification of the engine even-harmonic firing frequencies. This 
possibility permitted a muffler design with a less demanding transmission-loss perform- 
ance at the engine-bank fundamental firing frequency. This relaxed performance require- 
ment is evidenced by the 9 dB desired minimum transmission loss at the engine-bank 
fundamental firing frequency indicated by the arrow in figure 8. 

Back-pressure considerations . - Engine back pressure produced by the muffler sys- 
tem was not computed or measured; however, an effort was made to minimize back pres- 
sure by making the total exhaust-pipe area as large as was practical. Also, the inlet and 
outlet pipes in all the chambers were arranged so that the center lines coincided insofar 
as was practical. Finally, sudden flow-direction changes were avoided as much as 
possible. 
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Muffler configuration . - The computer -generated component lengths and a sketch of 
the corresponding geometrical configuration are shown in figure 13. It is important to 
understand that the general configuration is determined by the designer. The designer 
must specify constraints on the overall lengths, cross-sectional areas, relative arrange- 
ment and number of chambers, and the handling of the internal flow. The computer pro- 
gram adjusts the various component lengths within imposed constraints to achieve as 
nearly as possible a desired transmission-loss characteristic. The configuration shown 
in figure 13 produces the transmission-loss characteristic shown by the dashed line in 
figure 9, which compares favorably with the minimum desired transmission-loss char- 
acteristic shown by the solid line in figure 9. The muffler system was designed to pro- 
duce an exhaust-noise reduction of 15 to 20 dB(A). 

Installation .- The installation of the Y-connector and muffler is shown in figures 14 
and 15, respectively. The exhaust gases were combined with the Y-connector then led 
into the muffler inlet, as shown in figure 15. This particular mounting arrangement 
allowed minimum disturbance of the aircraft weight balance and prevented operating 
problems relative to the high-temperature exposed surfaces. The muffler volume was 
0.057 m^ (2.0 ft^) and the weight of the entire system was 21.3 kg (47 lb). 

Results and Discussion 

Exhaust-noise reduction .- Total radiated noise from the helicopter was measured 
with the muffler system installed (hereinafter called the modified helicopter) as was done 
for the standard helicopter configuration. The narrow-band (4-Hz) spectrum of the noise 
from the modified helicopter is shownin figure 16 by the solid curve, and the narrow-band 
spectrum for the standard helicopter from figure 7 is shown superimposed in figure 16 by 
the dashed line to facilitate evaluation of the muffler-system effectiveness. Clearly, the 
dominant components of the exhaust noise have been eliminated. At the low frequencies, 
rotor noise now appears to dominate the spectrum. In the higher frequency range, 320 to 
520 Hz, there is a 3 to 5 dB increase in the spectrum floor which may be attributable to 
self-noise (see ref. 16) generated by the muffler system. The overall exhaust-noise 
reduction obtained was approximately 11 dB(A), whereas the desired reduction was 15 to 
20 dB(A). 

Estimated muffler insertion losses .- The spectrum information shown in figure 7 
and similar data with the Y-connector only installed were used to estimate the insertion 
losses of the muffler at the engine-bank firing frequencies. These insertion losses are 
listed in the last column of table II. Also listed are the measured sound pressure levels 
under comparable test conditions corresponding to the standard exhaust system, 
Y-connector only, and Y-connector with muffler. From these tabulations the estimated 
muffler insertion losses were obtained. It is important to note that the listed insertion 
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losses are lower bounds on the actual muffler insertion losses because of possible con- 
tributions from other discrete-frequency noise sources radiating at approximately the 
same frequencies as the tailpipe exit. 

Estimated transmission losses .- If at a particular frequency the amplitude of the 
incident wave in the exhaust pipe outlet is not changed by the addition of the muffler, then 
the insertion loss is equal to the transmission loss at that frequency. Assuming this con- 
dition to hold, the estimated insertion losses of table II were equated to muffler transmis- 
sion losses and plotted on the computed transmission-loss curve for the muffler. The 
results for the first four engine-bank firing frequencies are shown in figure 17. For fre- 
quencies 2, 3, and 4 the computed transmission losses are high by 5 to 15 dB. For the 
higher engine-bank firing frequencies, the estimated insertion losses are negative, indi- 
cating that the muffler is amplifying at these frequencies by 9 to 14 dB. Plausible expla- 
nations for this behavior include both source -impedance effects and self-generated noise 
by the muffler at discrete frequencies (ref. 16). The mathematical model described in 
this paper fails to account for these effects. 

Engine performance loss .- Accurate back-pressure measurements for the modified 
exhaust system were not obtained; however, the aircraft operator did not observe any 
significant degradation of engine performance for hover flight conditions. 

CONCLUDING REMARKS 

An improved design method for expansion- chamber mufflers has been described 
and applied to a helicopter exhaust-noise problem. The resulting three-stage expansion- 
chamber muffler together with the exhaust-pipe Y-connector for combining the exhaust 
gases from all cylinders reduced the exhaust noise by 11 dB(A). Experimental compari- 
sons between the noise levels for the standard and modified exhaust systems indicated 
that the muffler was generating self-noise and/or interacting with the source impedance 
to result in an amplification of 9 to 14 dB at the fifth and sixth engine-bank firing 
frequencies. 

The muffler volume was 0.057 m3 (2.0 ft3) and the weight was 21.3 kg (47 lb), 
which were within the limits for a flightworthy system. There was no significant loss of 
engine performance as judged by the aircraft operator. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 12, 1973. 
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APPENDIX A 


COMPUTER PROGRAM EXRSIL 


Language : FORTRAN 

Purpose : To compute transmission loss at 151 values of wavelength for multiple -stage 
expansion-chamber mufflers. The program makes a specified number of attempts to 
optimize the transmission loss in accordance with a desired minimum transmission- 
loss characteristic. The effect of mean flow is included in the program. 

Use : Data cards are prepared to be read by statements shown below: 

READ 8, N0, C 
READ 10, SI, S2, S3 
READ 10, AMACH 
READ 10, XYZ(1,1),XYZ(1,2) 

READ 22, NUMBER 
READ 10, G, (1),X(1),H(1) 

8 FORMAT (I10,F10.3) 

10 FORMAT (3 F10. 5) 

22 FORMAT (15) 

Input variables are 

N0 Number of optimization attempts 

C Sound speed in exhaust gas of proposed muffler, ft/sec 

51 Tailpipe cross-sectional area, ft^ 

52 Annular cross-sectional area, ft^ 

53 Chamber cross-sectional area, ft^ 

AMACH Mach number of mean flow in tailpipe of muffler 

XYZ(I,1) = 1=1, 2, 3 Wavelengths, ft (minimum, intermediate, maximum) 

XYZ(1, 2)1=1, 2, 3 Desired minimum transmission losses at above wavelengths, 

dB 

NUMBER Number of stages in muffler system (should not exceed 4) 
G(I),X(I),H(I) 1=1, 2, 3, 4 Minimum, initial, and maximum lengths of muffler-stage 

components starting at tailpipe and working toward 
inlet, ft (note that there will be four lengths per 
chamber) 

Restrictions : Tailpipe Mach number should be in the approximate range 0.05 to 0.15. 

The ratio of expansion- chamber cross-sectional area to tailpipe cross-sectional area 
should not exceed approximately 6. The number of expansion -chamber stages should 
not exceed 4. 
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APPENDIX A - Concluded 


Method : The theoretical basis used in this program is given in references 3, 4, 5 and 
described briefly in this paper. 

Accuracy : Within the limitations of the theoretical model, the accuracy of the computed 
transmission losses is for practical purposes a function only of the accuracy of the 
input data. 

References : See references 3, 4, and 5 of this report. 

Run time : 200 to 600 sec 

Storage : 42,000s (CDC 6000 series) 
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END 
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TABLE I 


PHYSICAL CHARACTERISTICS OF THE STANDARD HELICOPTER 
Main rotor: 

Diameter, m (ft) 11.3 (37.13) 

Number of blades 2 

Blade chord, m (in.) 0.28 (11.0) 

Airfoil NACA 0015 

Blade area, m2 (ff2) 3.18 (34.27) 

Disk area, m2 (ft2) 100.8 (1085) 

Solidity 0.0314 

Tip speed, m/sec (ft/sec) 210.3 690 

Design operating speed, rpm 356 

Tail rotor: 

Diameter, m (ft) 1.7 (5.67) 

Number of blades 2 

Blade chord, m (in.) 0.10 (4.13) 

Blade area, m 2 (ft 2 ) 0.209 (2.25) 

Disk area, m 2 (ft 2 ) . . 2.35 (25.3) 

Design operating speed, rpm 1920 

General: 

Normal gross weight, kg (lb) 1111 (2450) 

Empty weight, kg (lb) 716.6 (1580) 

Overall length, m (ft) 13.29 (43.62) 

Power, kW (hp) 194 (260) 

Maximum level airspeed, knots 85 

Gear ratios: 

Main rotor mast 0.111 

Tail rotor 0.6 

Cooling fan 1.5 

Engine characteristics: 

Number of cylinders 6 

Operating cycle 4 

Exhaust-gas temperature at MAP of 48 cm Hg, K (°F) 1090 (1500) 

Gas-flow Mach number at MAP of 48 cm Hg and area of 0.0045 m 2 0.1 
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TABLE H.- ESTIMATED INSERTION LOSS OF MUFFLER OBTAINED FROM 
MEASURED SOUND PRESSURE LEVELS OF STANDARD AND 


MODIFIED EXHAUST SYSTEM 


Engine -bank 
firing-frequency 
harmonic 

Sound pressure level, dB, for - 

Estimated 
insertion loss 
of muffler, dB 

Standard 

exhaust 

system 

Y- connector 
only 

Y-connector 
and muffler 

1 

94.5 

81 

69.5 

11.5 

2 

83.0 

94 

69.5 

24.5 

3 

88.5 

65 

56.5 

8.5 

4 

77.0 

74 

64.2 

9.8 

5 

66.6 

54 

68.0 

-14 

6 

72.5 

60 

69.5 

-9.5 

7 

76.2 


69.5 


8 

70.0 


64.2 
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(a) Area expansion or contraction discontinuity. (b) Side branch discontinuity. 

Figure 2.- Simple duct discontinuities. 
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Figure 3.- Increase 
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Figure 4.- Schematic diagrams of area discontinuities. 
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Figure 6.- Comparison of transmission loss with and without flow. 
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Frequency, Hz 

Figure 7.- Narrow-band (4 Hz) spectrum of helicopter noise emission at 30.5 m (100 ft), an angle of 
90° to nose of helicopter, and a hover height of 1.83 m (6 ft). 
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Figure 12.- Test helicopter with standard exhaust system for one three -cylinder bank shown. 
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Figure 13.- Schematic of muffler configuration. 
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Figure 14.- Y-connector installation. 
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